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The effect of impurities on the isotope coefficient is studied theoreticaily in the framework of Abrikosov- 
Gor'kov approach generalized to account for both potential and spin- flip scattering in anisotropic super- 
conductors. An expression for the isotope coefficient as a function of the critical temperature is obtained 
for a superconductor with an arbitrary contribution of spin-flip processes to the total scattering rate and 
an arbitrary degree of anisotropy of the superconducting order parameter, ranging from isotropic s wave 
to d wave and including anisotropic s wave and mixed (s+d) wave as particular cases. It is found that both 
magnetic and nonmagnetic impurities enhance the isotope coefficient, the enhancement due to magnetic 
impurities being generally greater than that due to nonmagnetic impurities. From the analysis of the 
experimental results on LSCO high temperature superconductor, it is concluded that the symmetry of the 
pairing state in this system differs from a pure d wave. 



The isotope effect has played an important role in the development of phonon mediated mechanism of electron 
pairing in superconductors. Its discovery gave rise to the theory of phonon mediated mechanism of electron pairing 
by Bardeen, Cooper and Schrieffer (jj. The BCS theory gave the isotope coefficient a, defined by the relation 



T c cx M~ a , or equivalently, a = — fjf^f , equal to 1/2, in agreement with the experiments for mercury []2j. For simple 
superconducting metals like Hg, Zn, S, Pb, etc., the values of a were found to be very close to the BCS value 1/2. The 
deviations from the BCS theory found in superconducting transition metals and their compounds were reasonably 
well explained by taking into account the effects of Coulomb interactions ||], left out in the BCS theory, and by more 
realistic treatments based on Eliashberg equations Q] . 

The discovery of high temperature superconductors (HTSCs) brought a serious challenge to the BCS theory which 
had firmly established the phonon mediated mechanism of electron pairing as a dominant cause of superconductivity 
in most previously known superconductors. For conventional superconductors, both theory and experiment agree with 
the fact that a approaches the BCS value of 1/2 as T c becomes larger. However, this trend is violated in copper oxide 
HTSCs, where a has been found to vary with doping in different ways ||. For the optimum doping, corresponding to 
the highest critical temperature, the value of a is usually quite small compared to the BCS value of 1/2. For a certain 
doping level, e.g. in some La2Cu04 based HTSCs, the value of a becomes greater than 1/2 ||. Such a behavior of a 
has been considered as one of the strongest evidence for a new mechanism of superconductivity in HTSCs. However, 
this argument is not as strong as it appears. For example, it was shown that, within the BCS picture or within the 
more realistic Eliashberg approach, these deviations can be understood j(| by considering the effects of anharmonicity, 
energy dependence of the electronic density of states, pair breaking effects, isotope mass dependence of the carrier 
concentration, etc. 

The presence of impurities strongly affects various characteristics of HTSCs, including the isotope coefficient j?)-^). 
Earlier theoretical attempts to describe the isotope effect in impure HTSCs were based on either the Abrikosov-Gor'kov 
formula for T c of an isotropic s-wave superconductor containing magnetic impurities |^JlO| or the Abrikosov-Gor'kov- 
like formula for T c of an anisotropic superconductor that contains nonmagnetic impurities only [pr| , ^2[ . An increase 
in a with impurity concentration has been demonstrated, in qualitative agreement with the experiment. However, 
the theory predicted the universal dependence of the normalized isotope coefficient, a/ao, on the normalized critical 
temperature, T c /T c q, where ao and T c q are, respectively, the values of a and T c in the absence of impurities. This 
prediction does not agree with the experimental findings . It was shown by Kresin et al. |l(| that the universal 
behavior of a/ao versus T c /T c q is restricted to the case that magnetic impurities are the only cause for the increase of 
a/ao, and that other effects like the nonadiabaticity of the apex oxygen in YBa2Cu30y could break the universality. 
Here we provide an alternative explanation of how the nonuniversality of the dependence of a/ao on T c /T c q arises. 
Our approach is based on taking into account the combined effect of both nonmagnetic and magnetic scatterers on 
the critical temperature of a superconductor with anisotropic superconducting order parameter. 
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To derive the expression for the isotope coefficient, we make use of the equation for T c of an anisotropic supercon- 
ductor containing both nonmagnetic and magnetic impurities |L3]: 
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where is the digamma function; r^f is the electron relaxation time due to exchange scattering by magnetic impurities; 
r is the total electron relaxation time due to potential scattering by both magnetic and nonmagnetic impurities, as 
well as due to exchange scattering by magnetic impurities. Eq. ([[]) was obtained in Ref. |l3| within the weak-coupling 
limit of the BCS model in the framework of Abrikosov-Gor'kov approach. It generalizes the well-known expressions 
Ji4],[L5| for the critical temperature of impure superconductors to the case of combined effect of both nonmagnetic and 
magnetic impurities on the critical temperature of anisotropic superconductors. The coefficient x = 1 (^v)) fs 



quantifies the degree of anisotropy of the order parameter A(p) on the Fermi surface (FS), where the angular brackets 
(■■■)fs stand for a FS average. For isotropic s-wave pairing, (A(p))|, s = (A 2 (p)) F s, and hence x = 0- F° r a 
superconductor with d-wave pairing we have x = 1 since (A(p)) F s = 0. The range < x < 1 corresponds to 
anisotropic s-wave or mixed (d + s)-wave pairing. The higher the anisotropy of A(p) (e.g., the greater the partial 
weight of a d-wave in the case of mixed pairing), the closer to unity is the value of x- Note that in two particular cases 
of (i) magnetic scattering in an isotropic s-wave superconductor (x — 0) and nonmagnetic scattering only in a 
superconductor with arbitrary in-plane anisotropy of A(p) = 0, < x < I), Eq. ([j]) reduces to the well-known 

expressions [pr| , ^5| . 

It is convenient to specify the relative contribution of spin-flip scattering rate, l/r^f , to the total scattering rate, 
1/t, by the dimensionless parameter 7 defined as = 7^. The greater is the relative contribution from exchange 
scattering by magnetic impurities to 1/t, the higher is the value of 7 (7 ranges from in the absence of exchange 
scattering to 1 in the absence of non-spin- flip scattering). In general, the value of 7 depends on the scattering 
strengths of individual nonmagnetic and magnetic impurities, as well as on their concentrations. At relatively low 
doping level, one can expect 7 to depend only on the type of the host material and doping elements, not on the 
impurity concentration Jl^]. 

Differentiating Eq. (Ill) for T c with respect to the isotopic mass M under a reasonable assumption that electron 
relaxation times and the anisotropy coefficient x do not depend on M, taking the definition of the parameter 7 into 
account and using the definition of the isotope coefficient a, one has 
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Equation (g) is obviously more general than equations used previously for the analysis of the pair breaking effect 
on the isotope coefficient in HTSCs p|,|l0| |l2t . Indeed, Eq. (j|) does not rely on particular assumptions about the 
symmetry of the superconducting state and the nature of impurities (either magnetic or nonmagnetic). This equation 
can be used for an impure superconductor with arbitrary anisotropy of the superconducting order parameter and 
arbitrary relative contributions of potential and spin-flip scattering to the electron relaxation time. Such an approach 
is of particular importance for HTSCs doped with various chemical elements since, first, there is a strong evidence for 
a dominant d-wave (i.e., highly anisotropic) order parameter in HTSCs ]l6| with a subdominant s-wave component 
p7[ , and, second, a lot of experiments give evidence for the presence of magnetic scatterers (along with nonmagnetic 
ones) in doped HTSCs |§fl|. 

At given values of x and 7, Eqs. (|l]) and (En define the dependence of a/ao on T c /T c q. One can see from Eqs. (|l|) 
and (^|) that the dependence of a/ao on T c /T c q has a universal shape for both an isotropic s-wave superconductor 
(x = 0) with nonzero contribution of exchange scattering to the total scattering rate (0 < 7 < 1) and a d-wave 
superconductor (x = 1) with an arbitrary ratio of spin-flip and potential scattering rates (0 < 7 < 1), as it is seen 
from Fig. 1. 

Quite a different picture is realized in the case < x < 1> i- c -; f° r a mixed (d + s)-wave or an anisotropic s-wave 
superconductor. In this case the behaviour of a/ao as a function of T c /T c o essentially depends on the value of 7. As 
T c /T c q goes to zero, i.e., in dirty superconductors, the value of a/ao tends to 1/(1 — x) i n t ne absence of exchange 
scattering (7 = 0), while a/ao grows as a/ao ex T c q/T c in the case 7 7^ for all values of x- Thus, the universality 
of a/ao versus T c /T c q curve breaks down for < x < 1- This is consistent with experimental observations. Indeed, 
the studies of isotope effect in impure HTSCs show that the curves of a/ao versus T c /T c0 vary with the type of 
impurities, i.e., with the value of 7 (since different chemical elements contribute differently to spin- flip and potential 
scattering rates of charge carriers) . So, our results seem to be in a qualitative agreement with the experiment if one 
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assumes that the superconducting state in HTSCs differs from a pure s or d-wave, i.e., ^ / and ^ / 1, as is also 
indicated by several experimental observations fll7j , f22"| ] . 

However, our theoretical consideration doesn't account for a number of factors that could influence the isotope 
coefficient in impure superconductors, e.g., the energy dependent electronic density of states, nonadiabaticity, anhar- 
monicity, etc. In particular, the change in the carrier concentration upon chemical substitution can have a strong 
influence on T c (and hence on a) along with the pair-breaking effect. This is why, to compare our calculations with 
the experiment, we have taken the data on the isotope effect in the system Lai.ssSro.isCui-^MzC^ (M = Ni, Zn, Co, 
Fe), see Ref. M. In this system, the critical temperature decreases rapidly as the impurity content x increases, and 
falls down to « 0.3T c0 already at x = 0.02 — 0.03. At such low impurity concentration, the change in T c due to the 
change in can be neglected in the first approximation (as compared with the pair-breaking effect) since the value 
of T c q ~ 40K in the impurity-free HTSC Lai.ssSro.isCuC^ corresponds to the maximum on the curve T c (nh), and 
hence changes insignificantly with nh as long as the change in is small. Note that such an approach (i.e., ignoring 
the change in the carrier concentration upon doping) may not be appropriate in case of YBa2(Cui_ x Zn x )307 where 
T c /T c0 » 0.3 at x = 0.06 and all the more for Yi_ x Pr x Ba 2 Cu 3 7 where T c /T c0 w 0.3 at x = 0.5, see Ref. @. A 
reasonable explanation of the isotope effect in these HTSCs has been given by Kresin et al. jlo| who considered the 
interplay between the changes in the carrier concentration and nonadiabaticity of the apex oxygen. 

Fig. 2 shows the experimental data along with theoretical graphs calculated for x = 0-5 and different values of 7, 
ranging from to 1. One can see that at low impurity content x < 0.008 (T c /T c o > 0.75) there is a good agreement 
between the theory and the experiment. However, at higher values of a; (i.e., at lower values oIT c /T c q) the experimental 
points lie well above the theoretical curves for all impurity elements, except for M = Ni. The same is true for other 
values of \ since the upper curve in Fig. 2 (for 7=1) changes insignificantly with x, at least for T c /T c q > 0.2. 

To bring the theory closer to the agreement with the experiment, one can assume that the value of T c q in the 
impurity-free HTSC Lai.85Sr .i5CuO4 is strongly depressed relative to its "intrinsic" value Q because of the scat- 
tering of charge carriers by inhomogenities produced by substitution of Sr for La. Recently it was suggested |24| that 
the anomalous response of the anisotropic superconducting state to the development of low energy dynamical charge 
stripes can also cause the suppression of the "intrinsic" T c q. Taking T c q to be equal to its "intrinsic" value, it is 
possible to reach a qualitative agreement with the experiment on the isotope effect in impure HTSCs even for heavily 
doped samples with low T c . Since the value of "intrinsic" T c q in Lai.8sSro.i5Cu04 is not known a priori, we take the 
suggested value 24 of 90 K. 

Fig. 3 shows the results for the set of T c o = 90K and \ = 0.5. For such a choice of the values of T c q and %, the 
theoretical curves for 7 = — 1 are closer to the experimental data H . Since the value of ao at " intrinsic" T c q is not 
known, we assumed for simplicity that the value of a is the same at T c q w 40 K and 90 K. This assumption can be 
the reason for a discrepancy between the theory and the experiment in the region T c /T c q — 0.3 — 0.4, see Fig. 3, since 
one could expect the value of a at T c q = 90K to be somewhat lower than at T c q w 40K. It should be noted that in 
our theory different types of impurities correspond to different values of 7. The magnetic elements Fe and Co can be 
thought of being characterized by 7 in the range 0.1 - 1, while Ni, whose magnetic moment is reduced considerably by 
doping [pj, can be assigned somewhat lower value of 7 w 0.01, see Fig. 3. Similiarly, the doping by the nonmagnetic 
element Zn induces magnetic moments [|l8|,^5| . One can see from Fig. 3 that Zn can be assigned the value of 7 in the 
range 0.01 - 0.05, i.e. lower than in the case of Fe and Co but greater than in the case of Ni. 

Finally, we note that the agreement between the theory and the experiment becomes worse if one takes the value 
of x closer to unity. This implies that the symmetry of the pairing state in La-based HTSCs can differ considerably 
from a pure d-wave. As far as we know, up to now there were only indirect arguments in favor of rf-wave symmetry 
of the superconducting state in Lai.85Sro.i5Cu04 [E6[, while phase-sensitive experiments are unknown to us. The 
value of x ~ 0.5 is however close to that expected for the two-dimensional order parameter of the type A(k) = Ao 
[cos(/c x a)+cos(fcj,a)], see Ref. p7f . 

In conclusion, we have studied theoretically the effect of both magnetic and nonmagnetic impurities on the isotope 
coefficient in the framework of a generalized Abrikosov-Gorkov approach for the anisotropic superconductors. We 
have shown how the interplay between the potential and spin-flip impurity scattering gives rise to the nonuniversal 
dependence of a/ao versus T c /T c q in mixed (d + s) wave or anisotropic s wave superconductors. Our main result is 
that if the impurities are viewed as the only cause for the increase in the isotope coefficient in Lai.85Sro.i5Cui_ I M 2; 04, 
then the symmetry of the superconducting order parameter in Lai.ssSro 15CUO4 appears to be different from a pure 
ci-wave. Indeed, even if one takes into account relatively large experimental errors in the values of a, it is clearly 
seen that experimental points do not lie on a single "universal" curve of a/ao versus T c /T c q as one could expect 
in the case of a pure <i-wave symmetry of the superconducting state. According to our calculations, the difference 
in a/ao versus T c q/T c curves for different impurity elements can be attributed to different contributions from the 
exchange scattering to the total scattering rate of charge carriers in the mixed [d + s)-wave or anisotropic s-wave 
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superconducting state. The agreement with the experiment is much more better if one assumes that the "intrinsic" 
value of T c q in Lai.ssSro.isCuCU reaches « 90 K, more than twice greater than the experimental value w 40 K. It would 
be interesting to check if it is possible to explain the nonuniversality of a/ao versus T c /T c0 in La1.85Sro.15Cu1_zM.rO4 
within a pure d- wave framework. 
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FIGURE CAPTIONS 



Fig. 1. Universal dependence of the normalized isotope coefficient a/ao on the normalized critical temperature 
T c /T c q in an impure isotropic s-wave superconductor (x = 0) with a finite concentration of magnetic scatterers and 
an impure d-wave superconductor (x = 1) with an arbitrary ratio of spin- flip and potential scattering rates. 

Fig. 2. Same as in Fig. 1 for \ — 0.5 (a specific case of anisotropic pairing) for different values of the coefficient 7 
specifying the relative contribution to the total scattering rate from exchange scattering. 7 = (dot-dashed curve), 
0.01 (thin solid curve), 0.05 (dashed curve), 0.1 (dotted curve), 1 (thick solid curve). Experimental data from Ref. || 
for isotope effect in Lai.ssSro.isCui-^M^CU with different x and M = Ni (triangles), Zn (open squares); Co (closed 
circles); Fe (closed squares). Experimental values of T c as a function of x are normalized to the value of T c q = 37.5 
K at x = 0. 

Fig. 3. Same as in Fig. 2 for T c o = 90 K, see text for details. 
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